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We report the first observations of beam losses due to bound-free pair production at the interaction 
point of a heavy-ion collider. This process is expected to be a major luminosity limit for the Large 
Hadron Collider (LHC) when it operates with 208 p D 82 + i ons because the localized energy deposition 
by the lost ions may quench superconducting magnet coils. Measurements were performed at the 
BNL Relativistic Heavy Ion Collider (RHIC) during operation with 100 GeV/nucleon 63 Cu 29+ ions. 
At RHIC, the rate, energy and magnetic field are low enough so that magnet quenching is not an 
issue. The hadronic showers produced when the single-electron ions struck the RHIC beampipe 
were observed using an array of photodiodes. The measurement confirms the order of magnitude of 
the theoretical cross section previously calculated by others. 

PACS numbers: 29.20.Dh, 25.75.-q 



When fully-stripped heavy ions of atomic numbers Z\ , 
Zi are brought into collision at the interaction point (IP) 
of a collider, a number of electromagnetic interactions 
are induced by the intense fields generated by the coher- 
ent action of all the Z\^, charges in either nucleus (for a 
review, see Ref. [![). Some of these "ultra-peripheral" 
interactions have much higher cross sections than the 
hadronic nuclear interactions that are the main object 
of study. Among them, the bound-free pair production 
(BFPP), sometimes known as electron capture from pair 
production, occurs when the virtual photon exchanged 
by the ions converts into a pair, and the electron is cre- 
ated in an atomic shell of one of the ions: 



+ Z 2 + c^ 



(1) 



The resulting one-electron atoms have a slightly larger 
magnetic rigidity than the original bare nucleus. (Mag- 
netic rigidity is defined as p/(Qe) — Bp for a particle 
with momentum p and charge Qe that would have bend- 
ing radius p in a magnetic field B). Since the transverse 
recoil is very small, this "secondary beam" will emerge at 
a very small angle to the main beam. However it will be 
bent and focused less by the guiding magnetic elements 
and may be lost somewhere in the collider ring. 

It has long been known that this process, together with 
electromagnetic dissociation of the nuclei, could be a ma- 
jor contribution to the intensity and luminosity decay of 
ion colliders [3, Q • It was realized more recently [j, Q 
that, in certain conditions, the BFPP beam will be lost 
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in a well-defined spot, initiating hadronic showers in the 
vacuum envelope of the beam. The resulting localized 
heat deposition could induce quenches of superconduct- 
ing magnets. Detailed calculations have been given else- 
where for 2.76 TeV/nucleon 208 pb 82+ operation of the 
LHC at CERN @, i, 0, @. The consequent luminosity 
limit is expected to occur at a level close to the design 
performance. It is therefore vital to test our quantitative 
understanding of the features of the BFPP process in or- 
der to ensure safe operation of the LHC, uninterrupted 
by lengthy quench-recovery procedures. BFPP has been 
measured in fixed target experiments [t| [l(| EH at lower 
energy but not, until now, in a collider. An opportunity 
arose during 63 Cu 29+ operation of RHIC at Brookhaven 
National Laboratory. 

The flux of ions in the secondary BFPP beam emerging 
from the interaction point is the product of the collider 
luminosity and the BFPP cross section. The partial cross 
section for electron capture for a given low-lying atomic 
bound state i on nucleus 1 has the approximate form [l2j 



cr, « Z 1 5 Z 2 2 (^ l l0g7cm + 5. i 



(2) 



where 7 cm is the Lorentz boost of the ions in the centre- 
of-mass frame and A iy Bi depend only weakly on Z . The 
total cross section for BFPP to one of the colliding ions 
is then 



ctbfpp 



(3) 



where the sum is over all atomic shells. Some numerical 



values from 12| are given in Table fU The cross section 
for collisions of 100 GeV/nucleon 63 Cu 29+ was not cal- 
culated in Ref. \YM but we have estimated it as 0.2 barn 
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TABLE I: BFPP cross sections, typical peak luminosity, 
BFPP rates and relative changes in magnetic rigidity at RHIC 
and LHC. Values are taken directly where possible, or esti- 
mated by fitting sums of contributions of the form @, to 
the information in Ref. [l2T ]. 5 = 1/(Z — 1) is the fractional 
deviation of the magnetic rigidity. 
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by interpolation of the data given in Fig. 7 of Ref. [t2T ]. 
scaled with Z and log7 cm according to Eq. ((2|). Applied 
to other tabulated values in Ref. [12], this method pro- 
duces an agreement within 10%. 

The impact point of the lost ions is predicted by track- 
ing the orbit of the BFPP beam in the collider optics until 
it intersects the physical aperture (vacuum pipe) of the 
beam line. To lowest order in the uncorrelated betatron 
amplitudes, a, b, and relative magnetic rigidity deviation 
from the central value, 5, the horizontal displacement 
from the central orbit, at a distance s from the IP is 

x(s) = a^2/3(s) cos [/*(«)] + b \/2/3(s) sin ]/i(s)] + d x (s)S. 

( 4 ) 

Averaging over all ions gives the transverse emittance 
of the beam: (a 2 + b 2 ) = e; (3(s) is the usual "Twiss 
function" determined by the focusing properties of the 
accelerator, /i(s) = (3(u)~ 1 du is the betatron phase 
and d x (s) is the locally generated dispersion (d x (0) — 
d' x (0) = 0) that encodes the spectrometer effect of the 
bending magnets. 

An ion following the central trajectory (i.e. in the mid- 
dle of the bunch) in the BFPP beam enters the IP with 
a = b = S = and emerges, still with a = b = 0, but 
5 = 1/(Z — 1). This ion will hit the beam pipe at the 
first location where the horizontal aperture A x satisfies 
A x (s) = d x {s)/(Z - 1). In the arc, A x = 3.455 cm. 

During 197 Au 79+ operation at RHIC, 6 is not large 
enough for this to occur (x max ~ 2.5 cm). However, dur- 
ing 63 Cu 29+ runs [3, the larger S (Table [J) meant that 
the ions should be lost in a well-defined location where 
they can be detected directly. Thanks to the small cross 
section, there was no risk of quenching any superconduct- 
ing materials. Since the BFPP beam was not visible to 
the existing detectors, beam loss monitors in the form of 
PIN diodes (PDs) of the type Hamamatsu S3590 were 
mounted on the outside of the magnet cryostat down- 
stream of the IP of the PHENIX experiment. The PDs 
have a time resolution of 25 ns and the counts were read 
out every second. They were initially positioned in a 
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FIG. 1: (color online) Left: The energy deposition from a 
typical shower shown in a thin slice in the x — s— plane through 
the geometry. The red arrow indicates the impact point and 
the green and orange arrows show the positions of the PDs in 
the CPD and WPD. Right: The 3D model of the geometry 
as implemented in FLUKA around the impact point. The 
orange squares show the PDs in the WPD. 
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FIG. 2: (color online) The horizontal projection of the 1 a 
envelope of the BFPP beam and the nominal beam emerging 
from the PHENIX IP. The BFPP beam hits the inside of the 
vacuum chamber at 135.5 m. 

wide configuration (WPD), 3 m apart (see Fig. [1]), and 
later moved closer together (CPD), 0.5 m apart, around 
an observed count-rate maximum at 141.6 m. 

A linear interpolation of d x (s) between magnet ends, 
based on the design model of the collider, identified 
the elements within which the BFPP beam was lost: a 
9.7 m long dipole magnet starting at s = 129.6 m, fol- 
lowed by a 2.12 m drift space and a quadrupole mag- 
net, shown in Fig [T] More precise tracking, including 
chromatic effects ((^-dependences), and a final step of an- 
alytic orbit calculation inside the dipole, predicted an 
impact at s; ~ 135.5 m with an angle of incidence of 
(x'(si)) ~ 2.7 mrad. Tracking an assumed Gaussian dis- 
tribution of the amplitudes a and b in Eq. ^ gives the 
distribution of impact momenta in the spot around this 
point. Fig. [2] shows the main 63 Cu 29+ beam together 
with the 63 Cu 28+ BFPP beam propagating from the IP. 

To the extent that the RHIC hardware differs from 
the design model, the impact point may be shifted. The 
global distortion of the central orbit by misalignments of 
the magnets is partially corrected by localized corrector 
magnets. A least-squares fit was made to orbit data from 
beam position monitors (BPMs) and measured misalign- 
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ments (of order 1-2 mm) of the quadrupole magnets in 
the form 



'(«) = a; / 3(s) + ^G(s-s j ) 



x yj 0(s)/3(sj) sin \p(s) - ufa)] i 



(5) 



Here 0(s) is the unit step function, a;^ the first two terms 
in Eq. |[4]) and 9j are the angular kicks from misalign- 
ments and correctors at locations Sj. Unfortunately not 
enough data were available to make a satisfactory fit. 
Instead, the possible spread of the impact point was esti- 
mated by calculating the propagated error using Eq. {5} . 
A possible error on the horizontal quadrupole misalign- 
ments of 1 mm was assumed, corresponding to the 95 % 
confidence interval of the measurement and the fact that 
the magnets may have moved between early 2005 when 
the BFPP measurement was done and late 2006 when the 
misalignments were measured. The absolute error on the 
BPM data was taken to be 0.5 mm, reflecting the limited 
data and possible non-reproducibility. Combining these 
errors gives an impact point of 135.5 ± 2 m. Keeping the 
nominal optics, but displacing the beam pipe by 1.4 mm 
at the impact point according to the measured misalign- 
ment of the dipole, moves the impact point to 136.0 m. 

Several distributions of BFPP ions at the PHENIX 
IP, centered on different initial offsets and angles, were 
tracked to obtain their impact coordinates on the inside 
of the vacuum pipe. These were fed into a 3D model (see 
Fig. [T]) of the impact region geometry, the PDs and the 
dipole magnetic field, implemented in the FLUKA 2006.3 
Monte Carlo code 0, [l6| , in order to simulate the 
shower in the magnet and particles emerging from the 
cryostat. This code has been benchmarked by others in 
the relevant energy range flU ] . 

Within 2-2.5 m from the average impact point most of 
the shower inside the magnet has died out, although par- 
ticles outside propagate farther. This can be understood 
through a rough estimate: The typical average nuclear in- 
teraction length for the superconducting coil and the iron 
cold bore is about 20 cm. Convoluted with the impact 
point spot size of about 1.8 m, this gives a total shower 
length of about 2.8 m, assuming that the shower lasts 
five interaction lengths. The maximum PD signal was 
expected about 2.5 m downstream of the impact point. 
When the shower is contained within the dipole, moving 
the impact point along s only translates the shower. If 
the impact point is moved closer towards the end of the 
magnet or the angle of incidence decreased, more of the 
shower emerges into the void before the quadrupole and 
the profile changes qualitatively. A second peak in energy 
deposition appears at the entrance of the quadrupole and 
eventually exceeds the first one. 

The PDs count the number of minimum-ionizing par- 
ticles (MIPs) that pass their active area of 1 cm 2 . In the 
simulation, the number of MIPs entering each PD was 
recorded, assuming a detection efficiency of 30%. Since 
the PDs are much more sensitive to particles entering 
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FIG. 3: (color online) The PD signals from both configu- 
rations together with simulation results, averaged over the 
fills and normalized to a typical average luminosity of 9.1 x 
10 27 cm~ 2 s _1 . The simulation was produced from the central 
BFPP orbit in the nominal lattice but with the beam pipe 
misaligned by 1.4 mm at the impact location according to 
measurements. Due to large statistical error bars, the sim- 
ulation gives only the order of magnitude to expect in the 
measurement. 
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FIG. 4: (color online). Count rates measured on the Zero 
degree calorimeter (ZDC) luminosity monitors (black, right 
scale) and the three PDs with the highest signal (colors, left 
scale) during a store with the WPD. The data was binned in 
30 second intervals. A clear correlation between the luminos- 
ity and the PD count rates can be seen. 



through the side orthogonal to the cryostat, only these 
particles were counted. At the expected BFPP ion pro- 
duction rate of 4 kHz, and the 0.01 MIPs expected per 
PD per ion, pileup from multiple interactions was not 
an issue. Both WPD and CPD configurations were sim- 
ulated for various impact points and a typical result is 
shown together with measured data in Fig. [3] 

Data were collected from the PDs during 14 fills of 
RHIC with 63 Cu 29+ ions. An example of the recorded 
signal is shown in Fig. [4] together with the Zero Degree 
Calorimeter (ZDC) count rate (which is directly propor- 
tional to the luminosity). The PD signals show a good 
temporal correlation with luminosity and are very local- 
ized along s, close to the predicted impact position. With 
a 1 s binning, the correlation coefficient, r, between the 
PD signal and ZDC is 0.7 at 141.6 m and decreases grad- 
ually to 0.1 at PDs outside the predicted impact area. 
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FIG. 5: The ZDC count rate (proportional to luminosity, left) 
and background-subtracted PD count rate (right) vs. relative 

orbit displacement between beams in a van der Meer scan. A 

— t — ^ 2 

Gaussian fit of the form ^4exp( (x ^°> ) is shown together 
with the data. The two fits had a x — 0.49 mm and 0.44 mm 
and x 2 /DOF = 0.23 and 0.20, assuming a measurement error 
of 500 Hz on the ZDC signal and 0.3 Hz on the PD. 

With a 30 s binning the highest r rises to 0.98. This 
makes it very unlikely that some other source than the 
collisions could be responsible for the signals. 

This becomes even clearer from van der Meer scans in 
which the beam orbits are scanned transversely across 
each other and the luminosity variation is recorded as in 
Ref. Fig.[5]shows the ZDC rate and the signal of one 
of the PDs as a function of the orbit bump amplitude. 

Since the data were taken parasitically during normal 
colliding-beam operation, the PD signals were often pol- 
luted by other losses, e.g., from collimation. Therefore 
the cleanest data sets with least interference were picked 
out (near the beginning of each fill) and the averaged 
count rate with background subtracted, normalized to 
luminosity, was calculated for each PD (see Fig. [3]). The 
background noise level was calculated for each PD as the 
average signal without luminosity. 

In the WPD, the maximum signal was recorded by 
the PD at 141.6 m but moved to 140.5 m in the CPD. 
The count rates ranged from 1 to 20 Hz depending on 
PD position and luminosity. This agrees well with the 
simulation. However, the maximum in the simulation 



came from the PD at 138.6 m, which does not exactly 
agree with the measurements. This means that, in re- 
ality, the second peak in the energy distribution outside 
the cryostat is actually higher than the first, while in the 
simulation the first peak is the higher. As can be seen in 
Fig. [1] the peak in the shower could also escape between 
the PDs in the WPD. However, if the impact point is 
translated within the 2 m error bar, the simulated max- 
imum moves further downstream and a fair agreement 
can be found. Moreover, the PD signals themselves have 
error bars — apart from the statistical error a small rela- 
tive change in the counting efficiency between the PDs 
could change the result. This has not been measured so 
we give no numerical estimate. 

One might hope to extract a value for the cross 
section for BFPP between colliding copper ions at 
100 GeV/nucleon. However, because of the uncertain- 
ties in both the number of MIPs entering a specific diode 
in the simulation and also the recorded count rates, we 
can only conclude from the good agreement that the the- 
oretical estimate of the BFPP cross section has the right 
order of magnitude. 

No signal was detectable at 31 GeV/nucleon, consis- 
tent with the much lower (L x ctbfpp) (see Table U). 

In conclusion, we made the first measurements of the 
localized loss of a BFPP generated secondary beam in an 
ion collider. These measurements were done with copper 
beams at RHIC. We found the location of the maximum 
loss monitor signal at 140.5 m from the IP, within 2 m of 
its calculated location, and the measured event rates in 
our PD detectors within a factor of 2 of calculated ones. 
The deviations between calculated and measured values 
are consistent with estimated errors. This is a valuable 
test of our ability to make quantitative predictions of this 
effect for the LHC, where it is expected to be one of the 
most restrictive luminosity limits. 
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